Confronting the damping of the baryon acoustic oscillations with observation 
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We investigate the damping of the baryon acoustic oscillations in the matter power spectrum due 
to the quasinonlinear clustering and redshift-space distortions by confronting the models with the 
observations of the Sloan Digital Sky Survey luminous red galaxy sample. The chi-squared test 
suggests that the observed power spectrum is better matched by models with the damping of the 
baryon acoustic oscillations rather than the ones without the damping. 
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I. INTRODUCTION 



The baryon acoustic oscillations (BAO), the sound oscillations of the primeval baryon-photon fluid prior to the 
recombination epoch, left their signature in the matter power spectrum [l|, The BAO signature in the galaxy 
clustering has recently attracted remarkable attention as a powerful probe for exploring the nature of the dark energy 
component commonly believed to be responsible for the accelerated expansion of the Universe [3, 0, H, 0, Ij, 8]. The 
usefulness of the BAO to constrain the dark energy has been demonstrated 

I , and a lot of the BAO survey 

projects are in progress or planned [H, [H, [13, [Hj The BAO signature in the matter clustering plays a role 
of the standard ruler, because the characteristic scale of the BAO is well understood within the cosmological linear 
perturbation theory as long as the adiabatic initial density perturbation is assumed. 

However, the comparison of the BAO signature with observation is rather complicated. The observed galaxy power 
spectrum is contaminated by the nonlinear evolution of the density perturbations, the redshift-space distortions and 
the clustering bias. This enables us to use the galaxy power spectrum for other supplementary tests, in addition 
to the test of the expansion history of the Universe for the equation of state of the dark energy. For example, the 
redshift-space distortions probe the linear growth rate of the density fluctuations [13, EM 13 ■ The growth rate is now 
recognized to be very important as the test of gravity on the cosmological scales. 

In the paper pol | , some of the authors of the present paper investigated how the quasinonlinear density perturbations 
affect the BAO signature. Especially, we focused on the damping of the BAO signature. The semianalytic investigation 
on the basis of the third-order perturbation theory demonstrated that the BAO damping is sensitive to the growth 
factor Di{z) and the amplitude of the matter power spectrum a^. Here z is the redshift and the growth factor is 
normalized as Di(z) = a at a 1, where a is the scale factor normalized as a = 1 at the present epoch. As a result, 
a measurement of the BAO damping might be useful as an additional consistency test by enabling one to probe the 
growth factor multiplied by the amplitude of the matter perturbation Di(z)ag,. In the present paper, we extend the 
previous work to include the redshift-space distortions, and confront the BAO damping with the observed SDSS LRG 
galaxy power spectrum. Throughout this paper, we use units in which the velocity of light equals 1, and adopt the 
Hubble parameter Hq = 100/ikm/s/Mpc with h = 0.7. 



II. DAMPING OF THE BAO 



We start with reviewing the theoretical modeling of the BAO damping. The BAO signature is extracted from the 
matter power spectrum P{k, ii, z) at redshift z in the following manner, 

B{k,^J.,z)^^^^-l, (1) 

where ^ is the cosine of the angle between the line of sight direction and the wave number vector, and P(fc,/i, z) is 
the corresponding smooth spectrum without the BAO. As will be explained in detail below, we adopt the formalism 
developed by Matsubara [2l[ for theoretical modeling of P(fc,/i, z). The corresponding smooth spectrum P{k,fi,z) 
is computed in the same manner as P{k,fj,,z) but with the no- wiggle transfer function in Re f. [ij . As an alternative 
method, one can utilize the cubic spline fltting method to construct the smooth spectrum [22 |. which we adopt in 
comparison with observations. 
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The modeling of the quasinonlinear power spectrum has been investigated by many authors, based on both the 
perturbation theory and numerical simulations. As a nonperturbative approach beyond the standard perturbation 
theory, Matsubara proposed a model of the quasinonlinear matter power spectrum using the technique of resuming 
infinite series of higher order perturbations on the basis of the Lagrangian perturbation theory (LPT) [2l[. One of 
the advantages of using the LPT framework is the ability to calculate the quasinonlinear matter power spectrum in 
redshift space, which can be obtained by 

p(^pV(fc,M,^) = e-"('''^)^?(^)^W [Dl[z)Pl^{k,y,)+D\{z)P^'\k,^,)+a{^,,z)Dt[z)g{k)P^^[k,^^^^ , (2) 

where 

5(fc)=6^y^ d^Pl:li^)^ (3) 

P^^{k,^^) = {l + f^?fpl^^{k), (4) 

-^lin (^) linear matter power spectrum at the present epoch, / = dlnDi/dlna, and Q;(/i, z) = 1 + /(/ + 2)/i^. 

Also, P2^'* (^) is expressed as 

P^^'\k,^l) = pi^{k,^l) + P^^{k,^J), (5) 



where 



n.rn ^ ^ 1 \ ) 



P^l\k,f,)^il + ff,')Pl^{k)Y,f,'-r^ I drP^{kr)B„,n{r), (7) 







and Anmir, x) and i?„m(r) are given in Appendix B of Ref. [2l[. We take pj^pj,{k, ^, z) as P(fc, ^, z) in Eq. ([T]). 

Figure 1 shows the BAO signature. Except for the right lower panel, the dotted curve is the linear theory, while the 
solid curve is the result from the LPT formula at redshift z = 1 for = 0, 0.5, and 1, respectively, which is explicitly 
given by 

The right lower panel summarizes the n dependence, which is given by Eq. ([5]). Note that the case /x = is equivalent 
to the LPT formula in real space. The amplitude of the BAO signature is degraded compared with the linear 
perturbation theory. Thus the quasinonlinear clustering and the redshift-space distortions decrease the amplitude of 
the BAO. 

Let us introduce the correction function of the BAO damping W{k,fi, z) by 

Pap(fc, W{k, M, z)] Pii„(A:), (9) 

where i?iin(fe) is the BAO signature in linear theory. In the previous paper [20| . which was restricted to real space, it 
was demonstrated that the correction factor can be written in a rather simple form. One of the main results of the 
present paper is that a similar simple formula can be derived in redshift-space. After some computation similar to 
the one in [20| , we found that the leading factor of the correction function can be approximately written as 

mfc^,)-^_^i(f^^M(M no) 

^^^^^^'^ - l + a{,,z)Dl{z)m p^:^{K,) ^ 

where P22^{fc,/x) and P[;^^(fc,//) are defined as P22^{fc,/i) and P[;^^(fc,/i), respectively, but with the no- wiggle transfer 



function. The formula (|TOl) in the limit of /i = reduces to the previous result derived for real space [20| ■ The dashed 
curve in Fig. 1 shows the approximate formula ^ with PTI| . 

To demonstrate the validity of the approximate formula PTI)) . Fig. 2 shows the relative error |Pcx — 5ap|/|-Scx| at 
wave numbers PI, P2, P3, Tl, T2 and T3, which correspond to the peaks and troughs defined in Fig. 1, as a function 
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FIG. 1: The BAO signature according to the hnear theory and the LPT formahsm. Except for the lower right panel, the 
dotted curve is the linear theory, the solid curve is the LPT formula, and the dashed curve uses the approximate formula 
Eq. (|10p . for fi = 0, 0.5, and 1, respectively. The lower right panel summarizes the ^-dependence, obtained with Eq. ((Sjl. 
The (quasinonlinear) redshift-space distortion causes more damping of the BAO signature. Here the redshift is z — 1, and the 
cosmological parameters are h — 0.7, flm = 0.27, fit ~ 0.046, ris — 0.96 and as — 0.82 



of the redshift. The upper left panel is /i = 0, the upper right panel is /i = 0.5, and the lower left panel is /i = 1, 
respectively. The lower right panel is the result for the angular averaged power spectrum — ^ap°l/l^cx°l (^^^ 

below for details) . The approximate formula works at the 10 % level. 

Figure 3 shows the correction function W{k, fj,, z) as a function of the wave number k at redshift z = 1 for 
^ = 0, 0.5, and 1, respectively. It is obtained with the approximate formula pI7|) . The dotted-dashed curve is the 
correction function for the angular averaged power spectrum W^^°{k, z) (see below). Thus the BAO damping due to 
the redshift-space distortion is more efficient compared to the result in real space. 

In practice, the angular averaged power spectrum is used in measuring the BAO signature, which is expressed, as 
follows, using the power spectrum in the LPT formula: 

With the use of Eqs. ([T]) and we find that B^^°{k, z) is approximately written as 

i3,7(fc,z) - [1-W^^™(fc,z)]i?n„(fc), (12) 

with 

f^, duWik, u, z)P(k, u, z) 
W^^%k, z) = \ 1 . (13) 

d^P{k, /i, z) 
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FIG. 2; Relative error \SB/B\ = |_Bcx - -Bap|/|-Box| at the wave numbers PI, P2, P3, Tl, T2 and T3, which correspond to 
the peaks and troughs defined in Fig. 1, as a function of the redshift. Except for the right lower panel, the cases ^ = G, 0.5, 
and 1 are shown, respectively. The right lower panel shows the result for the angular averaged power spectrum \5B/B\ = 

I oave Dave I /I Dave I 
|-fc>ex - -Dap l/|-fc>ex |- 



The lower right panel of Fig. 2 shows the relative error IB^™ — B^\/\B^^°\ as a function of the redshift at the wave 
numbers PI, P2, P3, Tl, T2 and T3, which correspond to the peaks and troughs of the BAO. The dotted-dashed 
curve in Fig. 3 plots z) as a function of k at the redshift 1. 

Figure 4 compares the theoretical prediction of the LPT formula with the results from the iV-body simulations 
(30 reahzations). Each of our simulations used 512'^ particles in periodic cubes with side length 10'^ /i~^Mpc [23j . 
We apply a method to correct the deviation from the ideal case of infinite volume (see [1^ for details). The panels 
correspond to redshifts z — 3, 2, 1, and 0.5, respectively. One can see the agreement between the iV-body result and 
the theoretical prediction. 



III. COMPARISON WITH THE SDSS LRG POWER SPECTRUM 



Now we confront the theoretical predictions with observations. In particular, we use the SDSS LRG sample from data 
release 6. The SDSS data reduction procedure is the same as described in Ref. [^]. Here we utilize the cubic spline fit to 
consistently construct the smooth component for both the theoretical and the observational power spectra. However, 
the overall shape of the power spectrum from the LPT formula does not match the observational power spectrum. The 
power spectrum of the LPT formula shows the exponential suppression at large wave numbers, as is shown in Eq. ([2]). 
This feature can be understood as the nonlinear rcdshift-space distortion [2lj, the so-called finger-of-God effect, which 
is automatically formulated in the LPT formalism. This suppression factor matches a phenomenological model of the 
redshift-space power spectrum on large scales in Ref. [l^]. Such discrepancies probably arise from the truncation of 
higher order perturbations and ignoring the effect of galaxy clustering bias. This might make a systematic error in 
extracting the BAO consistently. To avoid this, we first construct the theoretical power spectrum by multiplying the 
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FIG. 3: The damping correction W{k,ii,z) at redshift z = 1 as a function of the wave number, for ^ = 0, 0.5, and 1, 
respectively. The dotted-dashed curve is the angular averaged case W'^^'^{k, z). Here the cosmological parameters are the same 
as those of Fig. 1. 



No. 
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0"8 


Us 


bo 


a[h-^Mpc^] 


xL 


^LPT 


-^simple 


-^lin — GOV 


XlPT-cov 


1 


0.27 


0.046 


0.82 


0.96 


1.7 


27.5 


15.1 


13.7 


13.7 


28.2 


23.0 


2 


0.27 


0.048 


0.82 


0.96 


1.7 


27.2 


13.9 


13.0 


12.9 


26.4 


22.1 


3 


0.27 


0.044 


0.82 


0.96 


1.7 


27.8 


16.6 


14.5 


14.6 


30.5 


24.1 


4 


0.27 


0.046 


0.80 


0.96 


1.75 


26.2 


15.1 


13.7 


13.7 


28.2 


23.1 


5 


0.27 


0.046 


0.84 


0.96 


1.65 


28.8 


15.1 


13.7 


13.7 


28.2 


22.9 
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0.27 


0.046 


0.82 


0.94 


1.7 


28.0 


14.9 


13.5 


13.5 


28.0 


22.7 
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0.27 


0.046 


0.82 


0.98 


1.7 


27.0 


15.2 


13.8 


13.8 


28.4 


23.2 
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0.26 


0.046 


0.82 


0.96 


1.65 


28.8 


15.3 


12.5 


13.0 


29.7 


22.3 


9 


0.28 


0.046 


0.82 


0.96 


1.75 


26.0 


15.3 


15.0 


14.5 


27.3 


23.8 



TABLE I: The results of the chi-squared test for the BAO signature for various cosmological models. Xlpt is based on the 
LPT power spectrum, while Xun assumes the linear power spectrum. Xsimpic the minimum chi-squared value in fitting the 
model B{k) = [1 — d»fc^]i3iin(fc). Xun-cov a-nd Xlpt-cov assume the linear and LPT power spectra, respectively, in evaluating 
Eq. (HSl). 



LPT power spectrum by the function fogc"*^ so as to match the SDSS LRG power spectrum, 

^fit(fc) = ble'^'^'pj^^^ik), (14) 

where bo and a are the fitting parameters. Figure 5 demonstrates the example Pfit{k), whose parameters are described 
in Table I, labeled as model no.l. One can see that this fitting function matches the observed power spectrum well. 
We also note that the BAO signature extracted using Pfit(fc) is not sensitive to the choice of bo and a. 

Figure 6 compares the BAO signatures extracted from the theoretical models and the SDSS LRG power spectrum 
of Fig. 5. We computed the chi-square as 

where B^^{ki) and B°^{k,) are the theoretical and observational BAO signatures at wave number ki, respectively, 
and AB{ki) is the error. In the computation, we used the data in the wave number range of 0.015 < k < 0.195. The 
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FIG. 4: Comparisons between the theoretical BAO signature and the resuhs from the A*'-body simulation (squares with error 
bars) a.t z — 3, 2, 1, and 0.5, respectively. The solid curve is the LPT formula and the dotted curve is the linear theory. 
The solid line is the BAO distorted by the nonlinear redshift-space distortions. The cosmological parameters adopted in this 
comparison are h = 0.7, Sim ~ 0.28, Sis = 0.046, Us = 0.96 and as = 0.82 



values of the chi-squared test for various cosmological models are listed in Table I. Here Xlpt result for the 

theoretical LPT model, while xfin that for the linear theory, which does not take the BAO damping into account. 
In this computation, we have not fitted any parameters, and the number of degrees of freedom is 19. Xlpt < xfin 
all of the models. This means that the models with the BAO damping match the observational results better. 

As an additional test, we compared the observational BAO signature with a very simple theoretical model B(k) — 
[1 — d1k'^]Biin{k), which includes the leading correction to the damping. Taking d■^, as a free parameter, we computed 
the chi-square and we found a minimum value at ~ Ah~^yLpc, for the models in Table I. Xsimpic minimum 
chi-square. Note that the case d^, — corresponds to the linear theory. Then, A^^ = Xun — Xsimpic ~ 1' which suggests 
that the detection of the BAO damping is at the 1 sigma level. 

To see the effect of the covariance between the data points, we compute 

xL = E [^''(^^) - ^°'(^')] Pih)Cov~\h,k,)Pik,) [B'\k,) - B°^ik,)] , (16) 

id 

where Cov{ki,kj) is the covariance matrix of the power spectrum. Here the covariance matrix is obtained by using 
100 mock catalogs generated via the second-order Lagrangian perturbation theory and Poisson sampling. The details 
of the procedure are described in Ref. Q . Figure 7 shows the resulting correlation matrix, which is defined by 

T{h, kj) = GOY[k,,k,) ^^^^ 

y/Cov{ki,ki)Cov{kj,kj) 

The result of Xcov value is shown in Table I, where Xlpt-cov^ (Xun-cov^) the result for the theoretical LPT model 
(the hnear theory). We find Xlpt-cov < Xun-cov ^'^^ ^^'^ models, again. 
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FIG. 5: Comparison of the theoretical power spectra and the SDSS LRG power spectrum. The dotted curve is the LPT formula, 
and the short-dashed curve is the linear theory. The solid curve is the power spectrum multiplied by the correction factor (I14|) . 
The long-dashed curve is the linear theory multiplied by the constant b^. The cosmological and fitting parameters are described 
in Table I labelled as model no.l. 




FIG. 6: Comparison of the BAG features extracted from theoretical models and from the observational data. The dotted curve 
is the linear theory and the solid curve is the LPT result. The squares with the error bars are the results from the SDSS LRG 
sample. The chi-squared score for this example is listed in Table I, labelled as model no.l. 
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FIG. 7; Correlation matrix I3| obtained from 100 mock catalogs. 

As discussed in Ref . [20j . to the leading order, the magnitude of the BAO damping is proportional to the amphtude 
of the matter power spectrum, Di{z)as- Thus, the precise measurement of the BAO damping might be useful in 
determining Di(z)ag,- To estimate the minimum achievable error we have computed the diagonal entry for Di{z)a^ 
in the inverse Fisher matrix. For the Fisher matrix calculation we have adopted the same approach as described 
in Ref. [20| but here used the angular average of P(fc,/i, z) instead of real space power spectrum. The results are 
almost the same as those for real space in Ref. [23|. The minimum attainable error of Di{z = 0.9)ct8 is ~ 0.1 x 

1/2 

(A^/2000deg.2) ' (at the 1 sigma level), where A A denotes the survey area. In this computation, we assume that the 
galaxy sample covers the redshift range 0.5 < z < 1.3, the mean number density of galaxies n — 5.0 x 10~^/i^Mpc^'^, 
and the clustering bias b = 2.0. Note that the error on Di{z — 0.9)ag, depends on the mean number density of galaxies 
and the clustering bias i2(| . 

IV. SUMMARY AND CONCLUSIONS 

In the present work we investigated the influence of the redshift-space distortions on the damping of the BAO in 
the matter power spectrum. The modeling was based on the work developed by Matsubara, which uses the technique 
of resuming infinite series of higher order perturbations within the framework of the Lagrangian perturbation theory 
[2l| . The result shows that additional BAO damping appears due to redshift-space distortions. We confronted the 
theoretical BAO signature with the observed power spectrum of the SDSS LRG sample. The chi-squared test suggests 
that the observed power spectrum favors models with the BAO damping over the ones without the damping. Though 
the statistical significance is not high, the BAO damping has likely been detected in the SDSS LRG power spectrum. 
In our modeling we have not taken into account the effect of the clustering bias on the BAO damping. This should be 
considered more carefully (cf. 25]); however, the authors of Ref. [1^ show that the BAO damping does not depend 
much on the halo bias in redshift space. 
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